Time evolution and decoherence of entangled states realized in coupled 

superconducting flux qubits 



in 
o 
o 

C 

a 

Hi 

oo 

(N 



13 
C 

o 

(N 
> 

OO 

in 

o 
in 
o 
-i— > 

c3 



c 

o 
o 



X 



Takuya Mouri 1,2 , Hayato Nakano 2,3 , and Hideaki Takayanagi 1,2,3 
1 Department of Physics, Tokyo University of Science, Tokyo 162-8601, Japan 
2 CREST, Japan Science Technology Agency, Saitama 332-0012, Japan 
3 NTT Basic Research Laboratories, NTT Corporation, Kanagawa 243-0198, Japan 

(Dated: February 2, 2008) 

We study theoretically how decoherence affects superposition states composed of entangled states 
in inductively coupled two superconducting flux-qubits. We discover that the quantum fluctuation of 
an observable in a coupled flux-qubit system plays a crucial role in decoherence when the expectation 
value of the observable is zero. This examplifies that decoherence can be also induced through a 
quantum mechanically higher-order effect. We also find that there exists a decoherence free subspace 
for the environment coupled via a charge degree of freedom of the qubit system. 

PACS numbers: 03.67.Lx, 74.50.+r, 85.25.Cp 



I. INTRODUCTION 

Josephson-junction circuits behave quantum mechani- 
cally, if they are sufficiently decoupled from their environ- 
ment. However, in reality, since a macroscopic quantum 
state is difficult to be decoupled from their environment 
completely, a quantum mechanical superposition state 
suffers decoherence. This is one of the central problems 
that must be solved before these circuits can be used for 
quantum information processing [1]. In order to minimize 
decoherence on quantum states, it is very important to 
search what is the dominant source of decoherence and 
how sensitively qubits feel the effect of environment. In 
this paper, we theoretically present that the quantum 
fluctuation of an observable in a coupled flux-qubit plays 
a crucial role in decoherence, even if the expectation value 
of the obseravable is always zero. We also report that we 
can extract a decoherence-free single-qubit basis from an 
inductively coupled flux-qubit system. 



II. INDUCTIVELY COUPLED TWO FLUX 
QUBIT SYSTEM 

Among qubits based on Josephson-junction circuits, a 
flux qubit was realized as a superconducting ring with 
three Josephson-junctions[2,3] and coherent oscillations 
were demonstrated [4, 5, 6]. We discuss decoherence ap- 
pearing in the dynamics of an inductively coupled two 
superconducting flux-qubits in Fig. 1(a) [7,8]. Four of 
Josephson-junctions have the Josephson energy Ej and 
capacitance C. The remaining two have aEj and evC, 
where a — 0.75 is a constant. The Hamiltonian of the 
system is given by 



H = 



C 



2e 



(ill + 7m + ih + im + a(i 2 L3 + im)) 



-£./(cos[7li] + cos[7biJ + cos[7 L2 ] + cos[7_r 2 ](1) 
+a(cos[7 L3 ] + cos[ 7jR3 ])) + ^ (A) j 2 M , 



where 7, [i = LI, Rl, L2, R2, L3, R3, M) is the phase dif- 
ference across each junction or the inductance. By calcu- 
lating the Hamiltonian with parameters Ej — 100GHz, 

E c = ^ = 8GHz, a = 0.75, M = ^mji^) 2 , we can 
obtain the energy dispersion shown in Fig. 1(b) and we 
call the point at / = 0.5 as the "degeneracy point". 



qubiH 



qubit2 




FIG. 1: (a)Effective circuit of inductively coupled two flux 
qubits. Two superconducting closed loops correspond to two 
flux qubits. A qubit has three Josephson junctions. An ex- 
tenal flux $1 or $2 is penetrating through each qubit loop. 
Qubits are coupled to each other through the mutual induc- 
tance M. (b)Energy dispersions obtained when identical ex- 
ternal fluxes are applied to the qubits, f = 4>i/$o = 
($0 = h/2e). 

We can approximate the Hamiltonian for the two-qubit 
system using Pauli matrices erf, erf, 

H s = hxal + A lC rf + h 2 (J z 2 



A 2 crf +jaf<T%. 



Here, subscripts 1 and 2 indicate the left and right 
qubits, respectively. The basis of the two-qubit state are 
I li-U),! Iit2>,| Til2> and I tits)) where | |i) corresponds 
to the state where the persistent current flows counter- 
clockwise and I J.,) to the clockwise state around the loop 
of the qubit i(i = 1,2). hi(i — 1,2) is the parameter 
that depends on the external flux $i for the qubit i, and 
hi = I p i{<S>i - = I P i<Z>o{fc - 5), where I pi {i = 1, 2) is 
the absolute value of the persistent current in the each 
qubit determined by Ej and a. $0 = h/2e is the flux 
quantum. Aj(z = 1,2) is the tunneling matrix element 

each qubit has, which is the order of y/EcEje~^ •> ' l E ° . 
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j is the coupling constant between two qubits, which is 
determined approximately by Ml v \I v i. If the external 
flux is equally applied {h\ = hi = h) and the two qubits 
are assumed to be identical (Ai = A 2 = A), we can 
simplify the Hamiltonian as 

H s = ho\ + Aof + ha z 2 + Aa z 2 + ja\al. 

By comparing energy eigenstates of the Hs and Fig. 
1(b) near the degeneracy point, we find that this approx- 
imated Hamiltonian Hg reproduces well the result of the 
original two flux-qubit Hamiltonian with j = 1.78GHz 
and A = 2.07GHz. Therefore, hereafter we analyze the 
decoherence in our coupled flux-qubit system with this 
Hs ■ In terms of the original flux-qubit system, erf corre- 
sponds to the flux the qubit i has, and of to the charge 
polarization in the junction of the qubit. At the degen- 
eracy point h = 0, we can obtain eigenvalues and eigen- 
states [7]. 

E = -vV+4A2, |^ o) = | Tlt2) + | 

-( W ^ +4A2 )(ltll2) + Ult2)), 

Ei = -j, = | Ti| 2 )-| I1T2), 

E 2 =j, IV2) = | TiT 2 )-| I1I2), 

E 3 = V.f+4 A2, |y> 3 ) = I T1T2) + I I1I2) (2) 
7 - Ji 2 + 4A 2 

- C 1 )(ITll2) + |llt2». 

Here, a remarkable fact is that all the eigenstates of 
the two-qubit Hamiltonian at the degeneracy point are 
strongly entangled ones. Adopting a state composed of 
these entangled states as an initial state, we will dis- 
cuss the time evolution of the density operator and some 
quantities. 



III. PHYSICAL MODEL OF DECOHERENCE 
AND QUANTUM MASTER EQUATION 



tion[9], 

1 f* 

Ps = J dti(i/(ti)[a;, [a;(-ti),ps]] 
-iv(h)[x,{x(-ti), p s }}), 

where 

v(t) = I dujcos[uut] coth[— — — ]J(lo), 
Jo 2fc B J 

poo 

r](t) = / dtv sm[Lot]J (lu) . 
Jo 

ks is the Boltzman constant and T is the absolute tem- 
perature. We adopt the spectral density J(u>) = g a "_j^ 2 ■ 
Here a is the cut off frequency and we set a = 5. 

IV. RESULTS AND DISCUSSION 

A. decoherence via flux degree 

First, we assume that the interaction between the sys- 
tem and the environment is caused only through the total 
flux of qubits; then we take x = a\ + erf . This is a rea- 
sonable situation because we can expect that the system 
suffers decoherence due to the coupling to the detector 
like a SQUID and the fluctuation of the external flux, 
both of which are related to small flux generated by the 
persistent currents flowing in qubits. Now, we consider 
a superposition state composed of the first and second 
excited states IV'12 (0)) = -j= (\tpi) + IV^)) at the degen- 
eracy point(/i = 0) as the initial state. We can analyti- 
cally calculate the time evolutions of expectation values 
(erf), (erf), (erf) and (erf) when there is no decoherence. 
Provided that it is possible to set independent probes 
which can detect both flux and charge on each of qubits, 
we can measure expectation values (erf), (erf), (af) and 
(erf) independently. |V>i2(i)) and the time evolutions of 
observables at h = are obtained as follows; 



In order to take into account decoherence, we con- 
sider the total system as the sum of the qubit system 
and the environment. The total Hamiltonian is ifrOT = 
Hs + He + -Hint, where Hs, He, and -Hint, are the 
Hamiltonian of the qubit system, the environment, and 
the interaction between them, respectively. We showed 
already Hs and the components of environment and in- 
teraction are described as follows; He = ^2 n Tiw n a\ l a n , 
Hint = x^2 k \k{au + a\.). We regard the environment 
as a boson bath, , a\ are the annihilation and the cre- 
ation operators of the environmental mode k. Afc is the 
strength of the coupling in each mode, x means a physical 
quantity of the qubit system that is coupled to the en- 
vironment. Assuming that the system-environment cou- 
pling is small enough, we can treat -Hint as a pertur- 
bative term so that we acquire a quantum master equa- 



|V>i2(*)) = cos[ji]| | x itx2) + isin[7'i]| Txilx2), (3) 

(iMt)kfliMt)) =0, 
(iMt)kf|iMt)> = o, 

(Vi2(t)K|Vi2(*)) =-cos[2jt], 

(Vl2(t)K|V>12(t)) =COs[2jt]. 

Here | \ xi ) = | U) + \ U), I Li) = -\ U) + I ii) 
(i = 1,2). This shows that this coherent oscillation is 
impossible to detect through the flux, but there is a 
possiblity to detect through the charge. We will show 
the result of the numerical calculation in the presence of 
decoherence, which is obtained by solving the quantum 
master equation with the above initial state. As we have 
expected in the analytical calculation, we realize that 
expectation values (erf) and (erf) obtained in numerical 
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calculations are always zero. The most interesting fea- 
ture of this example is that the expectation value of the 
physical quantity which interacts with the environment 
(erf+crf ) is zero. Although we would speculate that there 
is no coupling between qubits and the environment, this 
turns out not true according to the calculation results. 
Figure 2(a) shows the oscillations of (of) and (erf). We 







(a) 


1 1 151 A 1)10 

A A A/Wyv 


ry vy Vs v v 20 

V Time 




5 p4 10^ 



FIG. 2: (a)Time evolution of expectation values (erf) and 
(erf) in the presence of decoherence. (b)Time evolution of the 
eigenvalues pi(i — 1,2, 3, 4) of the density operator, (h = 0, 
j = 1.78, A = 2.07, g = 1) Time is normalized by A. The 
initial state is |-0i2 (0)) . Decoherence is caused through x = 

0\ + &2 ■ 

find that thier amplitudes are obviously suffering deco- 
herence. Time evolution of eigenvalues of the density 
operator in Fig. 2(b) helps us realize that a pure state 
is rapidly changing into a mixed state due to decoher- 
ence because the emergence of the finite second and later 
eigenvalues means the system becomes a mixed state. 
Therefore, we can make sure that, even if the expecta- 
tion value of the interacting quantity between system and 
environment is always zero, there exists certain effect of 
decoherence. 

Next, we discuss when the initial state is the super- 
position state IV'oi(O)) = ^TjdV'o) + l^i)) 01 t ne ground 
and the first excited states. Figure 3(a) shows oscilla- 
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FIG. 3: (a)Time evolution of expectation values (erf) and 
(erf) with decoherence via x = erf + erf. (b)Time evolution 
of eigenvalues of the density operator. The initial state is 

hMo)>. 

tions of the flux induced at the qubitl and qubit2, and 
these oscillations suggest that the possibility to detect 
the fluxes (erf), (erf) separately by setting an indepen- 
dent SQUlDs over each of the qubits. Figure 3(a) also 
shows that amplitudes of oscillations hardly decay and 
it seems that the oscillations are hardly affected by de- 
coherence. This suggests that it could be possible to 
realize coherent oscillation which is robust against the 



decoherence via x = a \ + erf . We can again observe the 
expectation value of the coupling between the qubit sys- 
tem and the environment is always zero ((af + <r|) =0). 
From the above two examples, surely we can make coher- 
ent oscillations where (Hint) = (x) J2k ^k(a>k + a\) = 
by using initial states which is composed by entangled 
energy eigenstates. However, one is strongly affected by 
decoherence and the other is not. Then, we are interested 
in considering how quantum fluctuations of the coupling 
x may become the source of decoherence [8]. In fact, the 
calculations obtaining Figs. 2 and 3 also give the finite 
quantum fluctuations 



V(M)\Sx^ 12 (t)) ~ 1.4, 
V(^oi(t)\Sx 2 \^ i(t)} ~ 0.75, 

where Sx 2 =x 2 — (x) 2 . It looks like that the quantum 
fluctuations of the coupling x and resulting magnitude of 
decoherence are determined by the curvatures of energy 
dispersions in Fig. 1(b) around the degeneracy point 
/ = 0.5. The quantum fluctuation of a; = a\ + erf 
qualitatively corresponds to the quantum fluctuation of 
an applied flux in the horizontal axis of Fig. 1(b). In 
fact, by a perturbation calculation under the condition 
\h\ <C |A|, \j\, the energy dispersion of an eigenstate \i) 
in the vicinity of the degeneracy point is given by 



Ei 



^-h 2 j2 



K +Crf|j0)(j0|crf + <Tf |*0> 



E? 



where E® and are the eigenenergy and eigenstate at 

the degeneracy point. So, the curvature ^pK* = 0, 1, 2) 
is approximately proportional to 

2<i0K +cr||j0)(i0|crf +a* 2 \i0), 
j 

that is the quantum fluctuation of a\ +<r|. As a result, 
we have reached a conclusion that the decoherence at the 
degeneracy point depends on the curvature of the energy 
dispersion around the degeneracy point. 

Now we discuss the concurrence which is defined as 
C{t) - V6" - V6 where &(i = 1,2,3,4) is the 
i th eigenvalue of A(t) = ^ i \pi\ 2 \{ipi\u\a\ ■ This is 
a very important measure because the concurrence tells 
us how the two qubits keep entanglement which provide 
us the possiblity to implement quantum parallelism cal- 
culations. Figure 4 shows the time evolution of the con- 
currence. When compared with the decays of the corre- 
sponding coherent oscillations in Figs. 2(a) and 4(a), we 
can find that the entanglement decays more rapidly than 
the observables (erf) and (erf )in the presence of decoher- 
ence. 



B. decoherence via charge degree 

Next, we suppose that the interaction between the sys- 
tem and the environment is only through the total charge 
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Concurrence 




FIG. 4: (a)Time evolution of the concurrence with the initial 
state 1^12(0)). (b)Time evolution of the concurrence with the 
initial state |^>oi(0)). 



freedom the qubits have; x = <rf + erf ■ For example, this 
type of interaction might appear in the presence of resid- 
ual charges near the qubit junctions. When the initial 
state is |^/>i2 (0)) = -^GV'i) + IV^)), as we have analyti- 
cally calculated, we have already known (af) and (erf) are 
always zero and (ipi2(t)\af\ipi 2 (t)) and (V>i2(*)|cf \ipi2(t)) 
oscillate sinusoidally. The numerical calculation of the 
time evolution gives the behavior of the expectation val- 
ues, shown in Fig. 5(a). We can see (af) and (erf) 
amazingly keep oscillating without any reduction of the 
amplitudes, and we have confirmed that the pure state 
is kept perfectly. We find this is because the two basis 
states and ^2) composing the initial state |V ; i2(0)) 
construct a Decoherence Free Subspace(DFS)[10] for the 
coupling between the two qubit system and the environ- 
ment via x = af + erf [10,11]. In other words, \ipi) and 
^2) are the degenerate eigenstates of x. Therefore, we 
can suggest that it is possible to design a single quantum 
bit which does not form entanglement with the environ- 
ment through the charge degree of freedom for a super- 
position state of and \tp2)- This is an example of a 
suppression of decoherence by introducing a redundant 
qubit and a qubit-qubit interaction, like quantum error 



correction codes[12]. Figure 5(b) represents time evolu- 
tion of the concurrence and we can make sure that the 
entanglement never deteriorates because of the suppre- 
sion of decoherence in the DFS. 
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FIG. 5: (a)Timc evolution of the expectation values (erf), 
(erf). (b)Time evolution of concurrence. The initial state is 
|V'i2(0)). Decoherence is caused through x — af + erf. 



V. SUMMARY AND CONCLUSION 

In summary, we analyzed the time evolution of an in- 
ductively coupled flux qubit system in the presence of 
environment. We confirmed the importance of quantum 
fluctuations of the observable coupled to the environ- 
ment. We also found two basis states constructing a DFS 
for charge fluctuations. 
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